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O ; ABSTRACT 

' Far-ultraviolet (FUV) photons expel electrons from interstellar dust grains and the excess 

kinetic energy of the electrons is converted into gas thermal energy through collisions. This 



O ' photoelectric heating is believed to be the main heating mechanism in cool HI clouds. The 

^ , heating rate cannot be directly measured, but it can be estimated through observations of the 

^ ' [CII] line emission, since this is believed to be the main coolant in regions where the photoelectric 

effect dominates the heating. Furthermore, the comparison of the [CII] emission with the far- 
infrared (FIR) emission allows to constrain the efhciency of the photoelectric heating, using model 
calculations that take into account the strength of the radiation field. Recent [CII] observations 
^ \ carried out with the ISO satellite have made this study possible. 

In this work we study the correlation between FUV absorption and FIR emission using 
three-dimensional models of the density distribution in HI clouds. The density distributions 
are obtained as the result of numerical simulations of compressible magneto-hydrodynamic tur- 
bulence, with rms sonic Mach numbers of the flow ranging from subsonic to highly supersonic, 
0.6 < Ms < 10. The FIR intensities are solved with detailed radiative transfer calculations. 
The [CII] line radiation is estimated assuming local thermodynamic equilibrium where the [CII] 
line cooling equals the FUV absorption multiplied by the unknown efficiency of the photoelectric 
heating, e. 

The average ratio between the predicted [CII] and FIR emissions is found to be remarkably 
constant between different models, implying that the derived values of e should not depend on 
the rms Mach number of the turbulence. The comparison of the models with the empirical data 
from translucent, high latitude clouds yields an estimate of the photoelectric heating efficiency of 
e '^2.9xl0~^, based on the dust model of Li & Draine. This value confirms previous theoretical 
predictions. 

The observed correlation between [CII] and FIR emission contains a large scatter, even within 
individual clouds. Our models show that most of the scatter can be understood as resulting from 
the highly fragmented density field in turbulent HI clouds. The scatter can be reproduced with 
density distributions from supersonic turbulence, while subsonic turbulence fails to generate the 
observed scatter. 1 
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— ISM: kinematics and dynamics 



1. Introduction 

In diffuse interstellar clouds {Ay < 1) the 
gas is heated mainly by photoelectrons expelled 
from dust grains by far-ultraviolet (FUV) pho- 
tons (de Jong 1977; Drainc 1978: Bakes & Tie- 
lens 1994; Weingartner & Draine 2001b) in the 
energy range 6eV< hv <13.6eV. The high energy 
limit corresponds to the cut-off in the FUV ra- 
diation field caused by the hydrogen absorption 
(/ii/=13.6eV), while the low energy limit corre- 
sponds to the energy needed to free electrons from 
the grains {hv ~6eV). In the cold neutral medium 
(2kin ^ 200 K) photoelectric heating accounts for 
most of the heating, the X-ray and cosmic ray 
heating rates being more than an order of magni- 
tude smaller (Wolfire et al. 1995). In a relatively 
dense neutral medium (uh ^100 cm~^), where a 
significant fraction of carbon is in the neutral form, 
carbon ionization becomes an important heating 
source, but it is still not comparable to the photo- 
electric effect. 

The efficiency of the photoelectric heating, e, is 
defined as the ratio between the kinetic energy of 
photoelectrons available to heat the gas and the 
FUV energy absorbed by the grains. Bakes & 
Tielens (1994) estimated the photoelectric yield 
to be close to 10%, meaning that of all the pho- 
tons absorbed by grains only one in ten ionizes 
them. Furthermore, since the average absorbed 
photon has an energy of ~8 eV and the work func- 
tion for a neutral grain is ~5.5eV, only one third 
of the absorbed energy goes into the kinetic energy 
of the electron. The overall efficiency is therefore 
expected to be close to 3%, at least in the cold 
neutral medium. The rate decreases as the grains 
become positively charged. In the intense FUV 
radiation field of photodissociation regions, for ex- 
ample, the efficiency can be more than a factor of 
ten lower than in the cold neutral medium. 

Ingalls et al. (2002) have recently studied a sam- 
ple of translucent, high latitude clouds (HLCs). 
The [C+] 158 /im line ^^3/2 P1/2 was observed 
with ISO and this was correlated with FIR inten- 
sities. In the cold neutral medium the [CII] line 
dominates the cooling (Wolfire et al. 1995), neutral 
carbon being the next most important coolant. In- 
galls et al. (2002) estimated that for these clouds 
the CI cooling rate is at most 30% of the [CII] cool- 
ing rate and the CO rate is more than one order 



of magnitude smaller. Therefore, for the observed 
HLCs the [0+] 158/xm line should account for at 
least 60% of the cooling and in most lines of sight 
the percentage should be considerably higher. 

To be a tracer of photoelectric heating the [CII] 
emission must originate from the same regions 
where the photoelectric effect is the main heating 
source and, conversely, all regions where the pho- 
toelectric heating is dominant must be cooled pre- 
dominantly by the [CII] emission. As discussed by 
Ingalls et al. (2002), both conditions are approx- 
imately satisfied. The ionization potential of car- 
bon, 11.3 eV, is above the work function of grains, 
and therefore C^ exists only in regions where pho- 
toelectric heating is possible. Conversely, in re- 
gions where C+ is not present and C and CO are 
significant coolants, the photoelectric heating rate 
must be much lower than in regions where C+ is 
present. As in Ingalls et al. (2002) we therefore as- 
sume that the above conditions are satisfied and 
the photoelectric heating is exactly balanced by 
the [CII] line cooling. 

Since all the FUV radiation in the energy range 
6eV< hv <13.6eV absorbed by grains contributes 
to the photoelectric heating, the efficiency e of this 
heating mechanism can be computed as the ra- 
tio of the observed [CII] line emission and the ab- 
sorbed FUV radiation. If dust properties are as- 
sumed to be known, the absorbed FUV radiation 
can be estimated through the observed FIR emis- 
sion. The problem of estimating e is therefore re- 
duced to the problem of computing the correlation 
between the absorbed FUV radiation and the ob- 
servable FIR emission. Such correlation depends 
on both th(^ intensity of the radiation field and 
the column density and can be determined with 
model calculations. The degeneracy between the 
effect of the cohimn density and of the intensity of 
the the external radiation field could in principle 
be broken using the information contained in the 
shape of the spectral energy distribution. This re- 
quires, however, quite precise modeling of the dust 
emission. 

Ingalls et al. (2002) modeled their observations 
using a plane parallel geometry. The FIR emission 
was calculated with a dust model consisting only 
of grains at equilibrium temperature. A compar- 
ison with [CII] and FIR observations resulted in 
an estimate of e = 0.043 for the efficiency of the 
photoelectric heating. The average radiation field. 
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X, was estimated to be above the local interstellar 
radiation field (ISRF), x ~ l-6xo, where xo refers 
to the sum of the 2.7 K cosmic background and 
the spectrum given by Mezger et al. (1982) and 
Mathis et al. (1983). 

In the present work we follow the same ap- 
proach as in Ingalls et al. (2002), in the sense 
of deriving theoretical relations between the FUV 
absorption and the FIR emission and then using 
the enipiric;al rc;lation between [CII] and FIR in- 
tensities to estimate the average efficiency of the 
photoelectric heating. However, our models of 
the FUV absorption and FIR emission is signif- 
icantly more detailed and realistic than in Ingalls 
et al. (2002). We compute the radiative trans- 
fer on cloud models based on three dimensional 
numerical simulations of compressible magneto- 
hydrodynamic turbulence. The density distribu- 
tion of these turbulent flows provides a realistic 
model for the density inhomogeneity of interstellar 
clouds (Padoan et al. 1997a,b, 1998, 1999; Padoan 
& Nordhuid 1999). The penetration of FUV radi- 
ation into the cloud depends on the cloud struc- 
ture. Inside an inhomogeneous cloud the intensity 
of short wavelength radiation is higher than in- 
side a homogenous cloud. Furthermore, density 
variations generate a significant scatter in the re- 
lation between the local FUV absorption and FIR 
emission. Our model can therefore be used to es- 
timate what fraction of the observed scatter can 
be attributed to the inhomogeneous nature of the 
density field. Other factors that could also con- 
tribute to the observed scatter, such as anisotropy 
in the radiation field or abundance variations, are 
not considered. Other improvements of our work, 
compared with Ingalls et al. (2002), are related to 
the employed dust model. We use the three com- 
ponent model of Li & Draine (2001). The FIR 
intensities are calculated using the method of Ju- 
vela et al. (2003) and the emission from transiently 
heated dust grains is included. 

2. The cloud models 

The density distributions of the models are the 
result of three dimensional simulations of super- 
Alfvcnic, compressible, magneto hydrodynamic 
(MHD) turbulence. Three models arc used, from 
three simulations with different values of the rms 
sonic Mach number of the flow, Ms = 0.6, 2.5 and 



10.0 in model A, B, and C, respectively. 

The simulations are carried out on a staggered 
grid of 250^ computational cells, with periodic 
boundary conditions. Turbulence is set up as an 
initial large scale random and solenoidal velocity 
field (generated in Fourier space with power only 
in the range of wavenumbers 1 < < 2) and 
maintained with an external large scale random 
and solenoidal force, correlated at the largest scale 
turn-over time. The initial density and magnetic 
field are uniform and the gas is assumed to be 
isothermal. Details about the numerical method 
are given in Padoan & Nordlund (1999). 

Experiments are run for approximately 10 dy- 
namical times in order to achieve a statistically 
relaxed state. The cloud models used in this work 
correspond to the final snapshot of each simula- 
tion. The value of Ms is varied in different ex- 
periments by varying the thermal energy. The 
initial rms Alfvenic Mach number remains un- 
changed. Ma = 10.0, from run to run. The 
volume-averaged magnetic field strength is con- 
stant in time because of the imposed flux conser- 
vation. The magnetic energy is instead amplified. 
The initial value of the ratio of average magnetic 
and dynamic pressures is (Pm)in/(^'d)in = 0.005 
for all runs, so all the runs are initially super- 
Alfvenic. The value of the same ratio at later times 
is larger, due to the magnetic energy amplification, 
but still significantly lower than unity (0.21, 0.14 
and 0.12 for models A, B and C respectively). The 
turbulence is therefore super Alfvenic at all times. 

Supersonic and super-Alfvenic turbulence of 
an isothermal gas generates a highly inhomoge- 
neous gas density distribution, with a density con- 
trast of several orders of magnitude. It has been 
shown to provide a good description of the dy- 
namics of molecular clouds and of their highly 
fragmented nature (Padoan et al. 1997a,b, 1998, 
1999; Padoan & Nordlund 1999). Translucent 
high latitude clouds (HLCs), such as the ones ob- 
served by Ingalls et al. (2002), have been stud- 
ied thanks to their CO emission (Magnani et al. 
1996). Despite their lower gas density and presum- 
ably larger fraction of gas in the atomic form, the 
clouds can still be approximately modeled with an 
isothermal equation of state. At densities above 
1 cm~^ the gas temperature is expected to remain 
mostly in the range 50-100 K (e.g. Wolfire et al. 
1995). The observed clouds are characterized by 
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supersonic turbulent motions like denser molecu- 
lar clouds (see e.g. Miesch et al. 1999). It seems 
therefore appropriate to use supersonic turbulence 
as a model not only of dense molecular clouds, but 
also of more diffuse clouds, such as the ones stud- 
ied in this work. A subsonic model (model A) is 
also used in this work for the purpose of study- 
ing the effect of varying Ms over a large range of 
values. 

To compute the radiative transfer the data 
cubes of the density distributions are scaled to 
physical units by fixing the length of the computa- 
tional grid, L = 6.0 pc, and the average hydrogen 
density, (nn) = 100 cm~^. However, the radiative 
transfer calculations apply to all values of L and 
(nn) satisfying the condition that their product is 
constant, L (tih) = 600 pc cm~^= 1.8 x 10^^ cm~^. 
The column density in the models varies between 
different lines of sight, because of the inhomoge- 
neous density field. Variations are particularly 
large in model C, where the range of column den- 
sity values is comparable to that found in the 
Ingalls et al. (2002) sample. This can be seen 
for example by comparing the FIR intensities in 
the models and in the observations (see Sect. 4). 
When results are compared with observations, the 
angular size of the models (or their distance) must 
be fixed as well, in order to take into account the 
spatial resolution of the observations. One cell, 
6.0/250 pc, is assumed to correspond to 0.25', set- 
ting the model clouds at the distance of 330 pc. 

3. Calculation of FUV absorption and FIR 
emission 

Radiative transfer calculations are needed to 

obtain estimates for the [CII] emission and the 
total FIR emission along different lines of sight 
through the model clouds. Calculations are car- 
ried out on the original grid of 250*^ cells. The 
external radiation field is assumed to be the local 
interstellar radiation field near the Sun (Mathis 
et al. 1983). The dust model is taken from Li & 
Draine (2001). This includes silicate grains (sizes 
a > 3.5 A), graphite grains (a > 50 A) and PAHs 
(from a = 3.5 A to ~ 50 A) with ionization state 
corresponding to physical conditions in the cold 
neutral medium (see also Juvela et al. 2003). 

Following Ingalls et al. (2002) we assume ex- 
plicitly a balance between the cooling by [CII] line 



emission and the photoelectric heating. The lat- 
ter is taken to be equal to the dust absorption in 
the energy range 6-13.6 eV multiplied by the effi- 
ciency of the photoelectric emission, e. The pene- 
tration of the external radiation field is simulated 
with Monte Carlo methods, and the absorption in 
the given energy range is stored in each cell. The 
emerging [CII] intensity is calculated by summing 
the corresponding emission on selected lines-of- 
sight. assuming the emission is optically thin. The 
efficiency of the photoelectric heating is unknown 
but since it is assumed to be constant troughout 
the cloud it can be taken into account later in the 
analysis. 

The dust emission was calculated using the 'li- 
brary method' discussed in Juvela et al. (2003). 
First a mapping between the strength of the local 
radiation field at a few reference frequencies and 
the resulting dust emission is established. The ra- 
diation field at the reference frequencies is then 
simulated and the emission is obtained using the 
existing mapping. The emission, including the 
contribution from transiently heated particles, can 
be estimated with an accuracy of a few per cent 
even in the case of large three-dimensional mod- 
els (see Juvela et al. 2003). The infrared emission 
from the clouds is estimated using this method. 
Ingalls et al. (2002) calculated FIR emission based 
on 60/xm and 100/^m values using the formula 



^^g J(100/xm)- 2.58/(60Mm) ^^^^_6 
MJysr-i 



erg 



cm^ s sr 

(1) 

The formula gives the FIR emission between 
42.5 /Ltm and 122.5 /itm i.e. the wavelength range 
covered by the IRAS 60/im and 100/im filters 
(Helou et al. 1985). In our models the formula 
was found to be accurate to ~ 2 — 3 per cent. 
However, the quoted FIR intensities are here cal- 
culated by direct integration over the computed 
spectra. 

Both the FUV absorption (equivalent to the 
[CII] emission divided by the efficiency, /cii/e) 
and the FIR intensity, /fir, are calculated to- 
wards three directions perpendicular to the faces 
of the cubic cloud. For each direction, maps of 
250 X 250 positions are generated. Results are com- 
pared with observations presented by Ingalls et al. 
(2002). The ISO LWS observations of the [C+] 
158/im line have a resolution of 71". The FIR data 
are based on IRAS ISSA maps with a resolution 
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of 4 — 5 '. Wc assume the cells of our models have 
an angular size of 0.25', and the computed [CII] 
and FIR maps are convolved to the resolution of 
the observations. 

4. Results 

4.1. FIR and [CII] emission in model 
clouds 

For each model the FUV absorption and the 
FIR intensity, /fir are calculated towards three 
directions perpendicular to the faces of the den- 
sity field data cubes. Since the FUV absorption 
is assumed to be equivalent to the [CII] line in- 
tensity, /cii, divided by the efficiency, e, in the 
following we will refer to the quantity icii/e (with 
e still unknown), instead of the computed FUV 
absorption. 

Fig. 1 shows the distribution of the ratio 

-^cii/(£^fir) in the three models. The figure is 
obtained at the original resolution (250x250 pix- 
els), and neither 7fir nor /cii/e are spatially con- 
volved. Model A is based on a flow with subsonic 
turbulence, with only small density fluctuations. 
This is reflected in the top panel of Fig. 1 show- 
ing a smooth distribution of Icu/ i^Ipm) in model 
A. The density contrast increases with increas- 
ing Ms. Model B, with Ms = 2.5, shows a more 
'clumpy' distribution of /cii/(£-^fir) (middle panel 
of Fig. 1). Model C, with Ms = 10, has clearly the 
most inhomogeneous and filamentary distribution 
(bottom panel of Fig. 1). 

In the case of the nearly homogeneous cloud 
A, both 7fir and /cii/e are determined mainly 
by the distance to the cloud surface. Since /cii/e 
is assumed to be equivalent to the FUV absorp- 
tion, /cii/e depends only on photons in the energy 
range 6-13.6 cV, while /fir depends on a much 
wider wavelength range responsible for heating the 
dust grains. The dust extinction increases towards 
shorter wavelengths, and this leads to a reduction 
in the ratio /cii/(e/FiR) inside the clouds. At the 
cloud surface the ratio /cii/(e/FiR) is above one, 
and it decreases to <0.8 inside the cloud, in all 
three models. The total range of values is wider 
in the more inhomogeneous clouds, with the low- 
est values reached in the highest density peaks. 
The average ratios are, however, very similar in 
all three models, (/cii/(e/FiR)) ~0.9. The aver- 
age value depends, of course, on the total column 






Aa [arcmm] 

Fig. 1. — The ratio /cii/ (e/fir) in the three model 
clouds A, B and C (from top to bottom). The 
maps correspond to observations towards the y- 
direction. No convolution has been applied to ei- 
ther /cii/e or /fir data. 
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Fig. 2. — The ratio /cii/ (c-^fir) as a function of the visual extinction along the line of sight. The three panels, 
a, b and c, correspond to the three models, A, B and C respectively. The figures show values calculated 
from the models towards the y-direction. The relations are shown for both the original data (dots) and the 
convolved data (crosses). For model A the indicated area has been plotted also on a larger scale. 



density and it would become lower in more opaque 
clouds. The models have similar column densities 
as the clouds in the sample of Ingalls ct al. (2002). 
and if the clouds are subjected to similar radiation 
field as the models the exact correspondence can 
be checked by comparing the FIR intensities. 

In Fig. 2 we plot the ratio /cii / (c-^fir) against 
the visual extinction of the corresponding line of 
sight. The average optical depth of the model 
clouds is below one, but the total range of Ay 
values depends on the degree of inhomogeneity. 
While in model A the visual extinction stays be- 
low the other models include a number of 
lines of sight with Ay in excess of 2™ (model B) 
and 3™ (model C). The figure shows the results 
for both unconvolved data and for the case where 
data have been convolved to the resolution of the 
observations (71" for [CII] and -4' for FIR data). 
The Ay data has not been convolved. This would 
correspond to a situation where the extinction is 
derived from photometry of an individual back- 
ground star, probing only a single line of sight in- 
side the larger beam employed in the other obser- 
vations. 

The negative correlation between extinction 
and the ratio /cii/(e/FiR) is clearly visible in 
model A (Fig. 2a). In model B the correlation 
can be seen in the unconvolved data, while in the 
convolved data the large difference between the 
resolution of the [CII] and FIR observations in- 



creases the scatter and the correlation becomes 
weak. At high extinction values the convolution 
increases systematically the /cii/(e-^FiR) values. 
The [CII] intensity depends on FUV absorption, 
and because of the high optical depth the [CII] 
emission is likely to satiirate near the column den- 
sity maxima. Since the intensity distribution is 
already flat the convolution has little effect on the 
peak [CII] intensities. The FIR emission is less 
saturated and, more importantly, the FIR inten- 
sities are convolved with a beam that is several 
times larger. As the the FIR peaks are smoothed 
away the ratio /cii/(e^FiR) increases. The effect 
is even higher in model C because of the stronger 
density contrast. However, in model C the Ay 
dependence is quite weak even in the unconvolved 
data. This can be attributed to other effects re- 
lated to the increased density contrast. For ex- 
ample, dense filaments and sheets are common, 
with the result that the Ay along one line of sight 
can be much higher than the average extinction 
by which the material is shielded towards other 
directions. In the figure this results in a larger 
scatter along the Ay axis. 

In a homogeneous model the extinction of the 
external radiation field proceeds in a predictable 
fashion, and the ratio between optical and FUV 
intensities only depends on the distance from the 
cloud surface. In an inhomogeneous cloud there 
is no longer a one-to-one relation between the 
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Fig. 3. — The convolved quantities /cii/e and /fir as a function of visual extinction. The three panels, a, b 
and c, correspond to the three models, A, B and C respectively. For model A the indicated area has been 
plotted also on a larger scale. 



shape of the spectral energy distribution and the 
intensity of the radiation. The FUV radiation 
penetrates deeper into the cloud than in the ho- 
mogeneous model, and the variation in the ra- 
tio between optical and FUV intensities becomes 
smaller. If the cloud consisted of optically very 
thick clumps, the incoming radiation would have 
exactly the same spectral energy distribution at 
the centre of the cloud as on the surface - only 
the intensity would be smaller. MHD models are 
somewhere between the two extremes, that is be- 
tween a homogeneous cloud and a cloud with opti- 
cally thick clumps. In different parts of the cloud 
the spectrum of the radiation field may be differ- 
ent even for the same 'effective extinction' and this 
increases the scatter along the Icii/{elFm)-sx.is. 

In Fig. 3 the dependence on the visual extinc- 
tion along the line of sight is shown separately for 
/cii/e and /fir- This shows explicitly the sat- 
uration of the [CII] emission, which at high col- 
umn densities is caused by the total absorption of 
incoming FUV intensity. The [CII] line itself is 
assumed to be optically thin. The effect is partic- 
ularly clear in model B. It is present also in model 
C, but it is partially masked by the scatter result- 
ing from the more inhomogeneous cloud structure 
and the effects of the convolution. 

Inhomogeneity is expected to lead to a decrease 
in total energy absorption (see Juvela et al. 2003) 
and, therefore, to lower the emitted intensities. 
The FIR emission in models B and C is respec- 



tively 2% and 5% lower than in model A. The 

FUV absorption, or the equivalent quantity /cii/e, 
is 3.6% and 9.0% lower in models B and C respec- 
tively than in model A. The larger variation of 
/cii/e relative to /fir is the result of a difference 
in the optical depths relevant to the two compo- 
nents. For optically thin radiation the absorption 
depends only on the total amount of absorbing 
material while for optically thick radiation the ab- 
sorbed energy is proportional to the area filling 
factor. Since the optical depth is higher for the 
FUV radiation than for the longer wavelengths 
heating the grains, the quantity /cii/e is more 
sensitive to the density contrast (filling factor) of 
the cloud. The equilibrium dust temperature for 
large grains, and consequently the FIR emission, 
depends on longer wavelengths that are partly op- 
tically thin. Therefore, in more clumpy clouds the 
reduction in /fir is more modest than in /cii/e- 
However, the value /cii/(e/FiR) decreases from 
model A to model C by less than 2%. The in- 
terpretation of observations is affected by several 
other uncertainties, such as the imprecise knowl- 
edge of the dust properties or the strength of the 
external radiation field (see Sect. 5). A variation 
of 2% is insignificant compared with these other 
uncertainties. 

4.2. Efficiency of the photoelectric heating 

In Fig. 4 the quantity /cii/e is plotted versus 
the FIR intensity, /fir, for the models A, B and C 
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Fig. 4. — /cii/e vs. /fir in maps calculated from the three model clouds towards the y-direction. The 
three panels, a, b and c, correspond to the three models, A, B and C respectively. The [CII] data has been 
convolved to a resolution of ^1.2' and the FIR to a resolution of ^ 4'. The slope of the least squares line 
going through the origin is given in each frame. For model A the data in the indicated area has been plotted 
separately on a larger scale. 



(left to right) observed along the y direction. We 
show the relations for both the original data and 
the data convolved to the resolution of the obser- 
vations discussed in Ingalls et al. (2002), that is 
~ 1.2' for [CII] and ~ 4' for the FIR observations. 
The relations are close to linear, and Fig. 4 shows 
the least squares lines fitted to the convolved data. 
The range of intensities and the scatter in the rela- 
tions both increase from model A to model C. The 
slope of the relation decreases, albeit very slowly 
with increasing Mach number of the model. The 
fitted slopes for each of the three orthogonal di- 
rections of the lines of sight are given in Table 1. 

In Sect 4.1 the [CII] intensities were found to 
saturate at high column densities when all incom- 
ing FUV radiation was absorbed. As a result, at 
high intensities the points for unconvolved data 
fall below the least squares line. However, at the 
location of intensity maxima the convolution de- 
creases the FIR intensities more than the [CII] in- 
tensities, and in the figure the convolved points are 
shifted to the left more than down. This is seen 
especially in model C that has the largest density 
contrast. In model C the relation /cii/e versus 
JpiR for the convolved data turns slightly upwards 
at higher intensities, relative to the relation for 
the unconvolved data. In model B the convolution 
raises the points at high- Ay only slightly closer to 
the fitted least squares line. 



A comparison between the ratios /cii/(e^FiR) 
obtained from the models and the empirical re- 
lation /cii//fir allows to estimate the efficiency 
of the photoelectric heating, e. In Ingalls et al. 
(2002) the average ratio of CII and FIR intensities 
is found to be /cii//fir=(2.5±0.9)x10-2, for a 
sample of high latitude clouds observed with ISO. 
The slope they derive for large scale correlation us- 
ing DIRBE and FIRAS observations for |6| > 5 deg 
is also very similar, /cii=(2.54±0.03) x 10~^7fir - 
(1.1365±0.0004)x10"'5 ergcm'^ s'^ sr"!. 

We have calculated the efficiencies of the pho- 
toelectric emission, e, assuming Icu/Ifir = 2.5 x 
10~^. The results are given in the fourth col- 
umn of Table 1. We obtain an average value of 
e = 2.88 X 10~^. There is very little variation be- 
tween the models or the direction from which the 
model clouds are observed. However, as seen in 
Fig. 1, the variation within individual maps can 
be quite significant. In the case of unconvolved 
data in Fig. 1 the variation ranges from ~40% in 
model A to almost 90% in model C. 

5. Discussion 

The observations presented by Ingalls et al. 
(2002) contain scatter in the plot of /cii versus 
/fir exceeding the uncertainty of individual ob- 
servations. Differences between clouds can be at- 
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Table 1 

The slope of the relation /cii/e versus /fir in the model clouds observed from three 

DIFFERENT DIRECTIONS. THE VALUES ARE BASED ON CONVOLVED DATA. THE EFFICIENCY OF THE 
PHOTOELECTRIC HEATING, £, IS OBTAINED ASSUMING THE EMPIRICAL RELATION /cii//fIR = 2.5 X 10^ 

FROM Ingalls ET AL. (2002). 



Model 


direction 


< /cil/(e/FIR) > 


100 X e 


A 


X 


0.876 


2.85 


A 


y 


0.873 


2.86 


A 


z 


0.872 


2.87 


B 


X 


0.869 


2.88 


B 


y 


0.861 


2.90 


B 


z 


0.870 


2.87 


C 


X 


0.862 


2.90 


C 


y 


0.862 


2.90 


C 


z 


0.858 


2.91 



tributed largely to differences in the dust compo- 
sition. The photoelectric heating and therefore 
[CII] emission depends on the abundance of very 
small dust grains and PAH molecules while the 
far-infrared emission is caused mainly by large 
dust grains. The ratio between the abundances 
of very small grains and large grains is known to 
show significant spatial variations. The variations 
are more noticeable in dense clouds (e.g. Laureijs 
et al. 1991; Juvela et al. 2002; Stepnik et al. 2003), 
but changes have been observed even in cirrus type 
clouds (Bernard et al. 1999; Cambresy et al. 2001). 
Variations in the radiation field are probably less 
important. 

In the first approximation an increase in the 

strength of the radiation field affects FIR and [CII] 
intensities in the same way. However, as the in- 
creased dust temperature moves dust emission to 
shorter wavelengths the FIR intensity is reduced 
relative to the [CII] emission (FUV absorption). 
Assuming that the FIR intensity in the models is 
fixed and constrained by the observations, if the 
interstellar radiation field is increased the dust col- 
umn density must be decreased to maintain the 
FIR intensity unchanged. As a result of the de- 
creased column density, the [CII] emission is less 
saturated, and the ratio Icu/Ifir is increased. 
Variations within individual clouds could similarly 
be caused by variations in the local radiation field 
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IpiR ^ SI'S ^ s ' sr 



Fig. 5. — Comparison of the [CII] and FIR inten- 
sities in model C and in the sample of HLCs in 
Ingalls et al. (2002). For model C the points cor- 
respond to convolved data, and to an efficiency of 
the photoelectric heating of 2.88x10"^ 
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or in the properties of the gas or dust. 

In the present models dust properties are as- 
sumed to be constant throughout the clouds. The 
fractional abundance of C"*" is also kept constant 
and the [CII] emission is assumed to be responsible 
for all of the gas cooling. Even with these simpli- 
fying assumptions the computed (/ciIj^fir) values 
contain a significant scatter. The scatter is small- 
est in the most homogeneous model A. Only in a 
completely homogeneous model the points would 
fall on one line as both [CII] and FIR emission 
would depend solely on the distance from the cloud 
surface. The scatter is largest in model C, with 
the largest value of Ms and density contrast. In 
this model, both the range of Iqii and /fir val- 
ues and the scatter in their relation arc similar to 
the observations presented by Ingalls et al. Such 
a comparison is relevant only if the observed po- 
sitions were selected randomly within the clouds. 
This is not always the case, but the scatter of the 
observed points should still be close to the true 
variation within the observed clouds. 

In Fig. 5 we compare our model C with the ob- 
servation presented by Ingalls et al. (2002). We 
have plotted in the same figure our model C (di- 
rection y) for which the [CII] intensities have been 
calculated using the value e = 2.9 x 10^^ derived 
in Sect. 4.2. The model covers a large fraction of 
the observed range of Icu and /fir values. Af- 
ter subtracting a linear fit, the remaining scatter 
in all observations is about twice the scatter in 
model C, even when the observational errors are 
taken into account. On the other hand, in indi- 
vidual clouds the variation is very similar to that 
found in the model. Our results show that most 
of the observed scatter may originate from the in- 
homogeneous density structure of the clouds, and 
the observed scatter is consistent with that pre- 
dicted with supersonic and super-Alfvenic MHD 
turbulence with sonic rms Mach number Ms > 2. 
In subsonic models both the scatter and the to- 
tal range of [CII] and FIR intensities is instead 
significantly smaller than in the observed clouds. 

Compared with model C the scatter is signif- 
icantly larger only in the cloud MBM-12, which 
also has the highest FIR intensities. MBM-12 
has an HI column density of ~l.lxlO^^ cm~^, i.e. 
some 40% below the average column density of the 
models. Timmermann et al. (1998) suggested the 
cloud is subjected to stronger than average radia- 



tion field, X '^1^2x0, and this could also explain 
the high FIR intensities. On the other hand, since 
the cloud has significant molecular line emission 
(Magnani et al. 1996) the total column density 
may not be very different from the model, and 
therefore the observed FIR intensities do not im- 
ply a strong radiation field by themselves. Fur- 
thermore, if one would accept the hypothesis of a 
strong radiation field, the low FIR colour temper- 
ature (see Ingalls ct al. 2002, Fig. 4) would require 
a very low abundance of small dust grains. 

The scatter in the observed [CII] and FIR in- 
tensities may depend on two effects not included 
in the models: i) CI cooling and ii) optical depth 
of the [CII] line. Although the models assume all 
the cooling is from C+, CI emission may be im- 
portant on some lines of sight, increasing the ob- 
served scatter. According to Ingalls et al. (2002) 
the maximum contribution to the cooling from CI 
is approximately 30%. The scatter due to the the 
presence of CI cooling on some lines of sight should 
therefore be well below 15%. Furthermore, the CI 
contribution can be significant only for high col- 
umn density, and could not explain the scatter at 
low FIR intensities. The models also assume that 
the [CII] line is optically thin. If the line became 
optically thick at some lines of sight, the main ef- 
fect would be to flatten the /cii vs. /fir relation 
at high column densities, while the scatter would 
not be significantly affected. In the models, the 
saturation of Icu is caused by the total absorption 
of the FUV radiation that is converted into [CII] 
emission, not by self-absorption of the [CII] line 
itself. Since the observations do not show strong 
[CII] saturation with increasing FIR intensity, the 
optical depth of the [CII] line cannot be very high. 

The derived value of the photoelectric heating 
efficiency, e, is close to theoretically derived val- 
ues (de Jong 1977; Bakes & Tielens 1994). Bakes 
& Tielens (1994) performed theoretical calcula- 
tions using a MRN type grain size distributions 
(n ~ a~3-5^ Mathis et al. 1977). They found that 
about half of the heating rate is due to grains 
with sizes below 15 A and the contribution from 
grains larger than 100 A is insignificant. The 
net efficiency was calculated as a function of the 
parameter Go(r/K)^/^(ne/cm~^)~^, where Go is 
the intensity of the far-UV radiation field relative 
to Habing's field (Habing 1968), T is gas tem- 
perature and rie the electron density. The es- 
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timatcd efficiency is constant e ~ 3 x 10 ^ for 
Go(T/K)i/2(ne/cm-3)-i < 10^. This prediction 
should apply to the cloud sample in Ingalls et al. 
(2002). Bakes & Ticlcns found, however, that the 
predicted heating rate could be doubled if the ex- 
ponent of the grain size distribution was decreased 
from -3.5 to -4.0. Since Ay and FIR intensity de- 
pend on large grains and the heating mainly on 
PAHs the derived efficiency e is model dependent. 

In this paper we have used the dust model of Li 
& Draine (2001). Weingartner & Draine (2001b) 
based their studies on the same dust model and 
considered further different size distributions dis- 
cussed in Weingartner & Draine (2001a). For cold 
neutral medium (T ~100K) and size distribu- 
tions consistent with the diffuse cloud extinction 
law, Ry ~3.1, their predictions for the photoelec- 
tric heating rate are close to the values given by 
Bakes & Tielens (1994). The results depend on 
the carbon abundance which is limited by extinc- 
tion measurements to values below ~ 6 x 10~^ car- 
bon atoms per H nucleus (Weingartner & Draine 
2001a,b; Li & Draine 2001). The models of Li & 
Draine (2001) and Weingartner & Draine (2001a) 
favoured carbon abundances close to this upper 
limit and the resulting rate of the photoelectric 
heating is some 25% above the value given by 
Bakes & Tielens (Weingartner & Draine 2001b). 

Habart et al. (2001) derived values of e=2- 
3% for the efficiency of the photoelectric heating 
across the cloud LI 721. Compared with the In- 
galls et al. sample, the cloud has higher visual 
extinction. Ay ~ 3™, and, due to the proximity 
of a B2 star, it is subjected to a stronger radia- 
tion field, X ^ 5xo- By adopting the dust model 
of Desert (1990), Habart et al. were able to de- 
rive efficiencies of the photoelectric emission sep- 
arately for PAH. very small grain and large grain 
components. The derived efficiencies were ~ 3% 
for PAHs, ~1% for very small grains and ~0.1% 
for large grains. Part of the observed variation 
was attributed to changes in the abundance of the 
dust components. The clouds in the Ingalls sample 
have lower visual extinction, and abundance vari- 
ations are expected to be correspondingly smaller. 

Another source of uncertainty is the balance 
between the column density and the strength of 
the radiation field. As discussed above, if a 
stronger radiation field or a smaller column den- 
sity were adopted in the models, the calculated 



ratio /cii / (e /fir) would increase and the derived 
value of e decrease. According to Fig. 1, the ratio 
Icu/ (e/fir) increases up to 30% above its average 
value when moving from the cloud center to the 
cloud edges, where the radiation field is not atten- 
uated. Therefore, if the clouds had significantly 
lower extinction, the efficiency of the photoelec- 
tric heating could be lower by the same amount. 
We modified the model C by reducing its column 
density by a factor of three and by multiplying the 
intensity of the external radiation field by a fac- 
tor of two. The estimated average efficiency of the 
photoelectric heating decreased from 2.88xl0~^ 
to 2.25x10-2, i.e. only by 22%. Based on HI and 
molecular line data the column densities are in 
most of the observed clouds within a factor of two 
from the average column density of the models. In 
this paper we have used for the intensity of the lo- 
cal ISRF values given by Mezger et al. (1982) and 
Mathis et al. (1983). According to Draine (1978) 
and Parravano et al. (2003) the actual ISRF could 
be up to '--^70% stronger. However, the difference 
is very unlikely to be a factor of two. The actual 
errors caused by the uncertainty of the column 
density and intensity values is therefore less than 
22%. 

Compared with Ingalls et al. (2002) our es- 
timate of the photoelectric heating efficiency is 
smaller by one third. Most of the difference is 
due to a difference in the predicted FIR inten- 
sity. For example, for a model with optical depth 
Ay=1.0 mag through the cloud Ingalls et al. ob- 
tained a 100 /im surface brightness of slightly more 
than 15MJysr-^. For a homogeneous, spherically 
symmetric cloud with equal optical depth we ob- 
tain a surface brightness of ~9.5MJysr-^. This 
is close to what Bernard ct al. (1992) obtained 
for a similar cloud model using the dust model of 
Desert (1990). The difference between our result 
and Ingalls et al. cannot be due to model geome- 
try. One would expect more emission (higher /fir) 
from a three dimensional cloud heated from all di- 
rections than from a plane parallel cloud heated 
on two surfaces only. The most likely explanation 
for the high FIR intensities obtained by Ingalls et 
al. (2002) is their assumption of thermal equilib- 
rium. This means that all grains are constantly 
at temperatures close to 20 K and almost all ab- 
sorbed energy is re-radiated in the FIR. In our case 
a larger fraction of the absorbed energy is radiated 
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at shorter wavelengths because of the smaU par- 
ticles that are temporarily heated to much higher 
temperatures. 

6. Conclusions 

We have calculated FIR emission and [CII] 
line emission for three-dimensional density distri- 
butions of compressible magneto hydrodynamic 
turbulent flows with rms sonic Mach numbers 
Ms = 0.6, 2.5, and 10.0. The dust emission, /fir, 
is computed with full radiative transfer calcula- 
tions. The [CII] emission is estimated assuming 
the photoelectric heating caused by FUV photons 
between 0.0912/;im and 0.2066^m is balanced by 
[CII] line emission. The FUV absorption is deter- 
mined by the radiative transfer simulations, and 
is assumed to be equal to the [CII] line intensity 
divided by the unknown efficiency of the photo- 
electric heating, /cii/e- 

The average ratio /cii/(e/FiR) is in all models 
between 0.85 and 0.88, showing that its depen- 
dence on the density contrast (Mach number) is 
weak. However, dense filaments are visible in the 
maps as regions with lower value of /cii/C^/fir)- 
The degree of correlation between /cii/(e/FiR) and 
visual extinction decreases in more inhomogeneous 
clouds. In the case of simulated observations con- 
volved to different resolutions (71" for [CII] and 
~4' for FIR) most of the correlation is lost. 

The scatter in the observational (/FiR./cii/e) 
plot can be reproduced by models with rms Mach 
number Ms >2 (supersonic turbulence), showing 
that most of the scatter may be due to the inho- 
mogeneous nature of the clouds (likely due to the 
turbulence). In subsonic models (Ms <1) both the 
scatter and the total range of FIR and [CII] inten- 
sities become smaller than in the observed clouds. 
Using the empirical value 7cii//fir=2.5x10~^ 
found for high latitude clouds (Ingalls ct al. 2002) 
the efficiency of the photoelectric heating is found 
to be e ~ 2.9 x 10~^. The value is very close 
to the theoretical predictions for the cold neutral 
medium. 
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